Non-scattering wavenumbers and far field invisibility for a finite set of incident/scattering directions

We investigate a time harmonic acoustic scattering problem by a penetrable inclusion with compact support embedded in the free space. We consider cases where an observer can produce incident plane waves and measure the far field pattern of the resulting scattered field only in a finite set of directions. In this context, we say that a wavenumber is a non-scattering wavenumber if the associated relative scattering matrix has a non trivial kernel. Under certain assumptions on the physical coefficients of the inclusion, we show that the non-scattering wavenumbers form a (possibly empty) discrete set. Then, in a second step, for a given real wavenumber and a given domain D, we present a constructive technique to prove that there exist inclusions supported in D for which the corresponding relative scattering matrix is null. These inclusions have the important property to be impossible to detect from far field measurements. The approach leads to a numerical algorithm which is described at the end of the paper and which allows to provide examples of (approximated) invisible inclusions.Comment: 20 pages, 7 figure

Rayonnement sonore dans un écoulement subsonique complexe en régime harmonique (analyse et simulation numérique du couplage entre les phénomènes acoustiques et hydrodynamiques)

La thèse porte sur la simulation, en régime fréquentiel, du rayonnement acoustique en écoulement subsonique quelconque et dans un domaine infini. L'approche choisie s'appuie sur la résolution d'un système équivalent aux équations d'Euler linéarisées : le modèle de Galbrun. Ce modèle repose sur une représentation mixte Lagrange-Euler et aboutit à une équation dont l'unique inconnue est la perturbation du déplacement Lagrangien. Une des difficultés de l'approche de Galbrun est qu'une discrétisation directe de cette équation par une méthode d'éléments finis standard n'est pas stable. Un moyen de contourner cet obstacle est d'écrire une équation augmentée en ajoutant une nouvelle inconnue, le rotationnel du déplacement, appelée par abus vorticité. Cette approche conduit à un système qui couple une équation de type équation des ondes avec une équation de transport en régime fréquentiel. Et elle permet l'utilisation de couches parfaitement adaptées (PML) pour borner le domaine de calcul. La première partie du manuscrit est dédiée à l étude de l équation de transport harmonique et de sa résolution numérique, en particulier par un schéma de type Galerkin discontinu. Un des points délicats est lié au caractère oscillant des solutions de l'équation. Une fois cette étape franchie, la résolution du problème de propagation acoustique a été abordée. Une approximation basée sur l'utilisation d'éléments finis mixtes continus-discontinus avec couches parfaitement adaptées (PML) a été étudiée. En particulier, les caractères bien posés des problèmes continu et discret ainsi que la convergence du schéma numérique ont été démontrés sous certaines conditions sur l'écoulement porteur. Enfin, une mise en œuvre a été effectuée. Les résultats montrent la validité de cette approche mais aussi sa pertinence dans le cas d'écoulements complexes, voire d'écoulements dits instablesThis thesis deals with the numerical simulation of time harmonic acoustic propagation in an arbitrary mean flow in an unbounded domain. Our approach is based on an equation equivalent to the linearized Euler equations called the Galbrun equation. It is derived from a mixed Eulerian-Lagrangian formulation and results in a single equation whose only unknown is the perturbation of the Lagrangian displacement. A direct solution using finite elements is unstable but this difficulty can be overcome by using an augmented equation which is constructed by adding a new unknown, the vorticity, defined as the curl of the displacement. This leads to a set of equations coupling a wave like equation with a time harmonic transport equation which allows the use of perfectly matched layers (PML) at artificial boundaries to bound the computational domain. The first part of the thesis is a study of the time harmonic transport equation and its approximation by means of a discontinuous Galerkin scheme, the difficulties coming from the oscillating behaviour of its solutions. Once these difficulties have been overcome, it is possible to deal with the resolution of the acoustic propagation problem. The approximation method is based on a mixed continuous-Galerkin and discontinuous-Galerkin finite element scheme. The well-posedness of both the continuous and discrete problems is established and the convergence of the approximation under some mean flow conditions is proved. Finally a numerical implementation is achieved and numerical results are given which confirm the validity of the method and also show that it is relevant in complex cases, even for unstable flowsTOULOUSE-INSA-Bib. electronique (315559905) / SudocSudocFranceF

An alternative to Dirichlet-to-Neumann maps for waveguides

International audienceWe are interested by the treatment of the radiation condition at infinity for the numerical solution of a problem set in an unbounded waveguide. We propose an alternative to the classical approach involving a modal expression of Dirichlet-to-Neumann (DtN) operators. This method is particularly simple to implement since it only requires the solution of boundary value problems with local boundary conditions. The corresponding approximate solution is comparable in accuracy to the one obtained by truncating the infinite series in the DtN maps

Mathematical analysis of guided water waves

Projet IDENTIn this article, we present a detailed mathematical analysis of the phenomenon of water waves guided by the coast. We work with the mathematical model issued from the linearized theory of gravity waves. Our results mainly concern the existence, the number and the high and low frequency behaviours of the guided modes. These results point out a large variety of phenomena which are deeply influenced by the geometry of the coast. We obtain them with the help of the spectral theory of unbounded selfadjoint operators. Some numerical results are presented at the end of the paper to illustrate the theory

On the use of Perfectly Matched Layers at corners for scattering problems with sign-changing coefficients

International audienceWe investigate in a 2D setting the scattering of time-harmonic electromagnetic waves by a plasmonic device, represented as a non dissipative bounded and penetrable obstacle with a negative permittivity. Using the $\texttt{T}$-coercivity approach, we first prove that the problem is well-posed in the classical framework $H^1_{loc}$ if the negative permittivity does not lie in some critical interval whose definition depends on the shape of the device. When the latter has corners, for values inside the critical interval, unusual strong singularities for the electromagnetic field can appear. In that case, well-posedness is obtained by imposing a radiation condition at the corners to select the outgoing black-hole plasmonic wave, that is the one which carries energy towards the corners. A simple and systematic criterion is given to define what is the outgoing solution. Finally, we propose an original numerical method based on the use of Perfectly Matched Layers at the corners. We emphasize that it is necessary to design an $\textit{ad hoc}$ technique because the field is too singular to be captured with standard finite element methods

Iterative methods for scattering problems in isotropic or anisotropic elastic waveguides

International audienceWe consider the time-harmonic problem of the diffraction of an incident propagative mode by a localized defect, in an infinite straight isotropic elastic waveguide. We propose several iterative algorithms to compute an approximate solution of the problem, using a classical finite element discretization in a small area around the perturbation, and a modal expansion in unbounded straight parts of the guide. Each algorithm can be related to a so-called domain decomposition method, with or without an overlap between the domains. Specific transmission conditions are used, so that only the sparse finite element matrix has to be inverted, the modal expansion being obtained by a simple projection, using the Fraser bi-orthogonality relation. The benefit of using an overlap between the finite element domain and the modal domain is emphasized, in particular for the extension to the anisotropic case. The transparency of these new boundary conditions is checked for two- and three-dimensional anisotropic waveguides. Finally, in the isotropic case, numerical validation for two- and three-dimensional waveguides illustrates the efficiency of the new approach, compared to other existing methods, in terms of number of iterations and CPU time

A complex-scaled boundary integral equation for time-harmonic water waves

This paper presents a novel boundary integral equation (BIE) formulation for the two-dimensional time-harmonic water-waves problem. It utilizes a complex-scaled Laplace's free-space Green's function, resulting in a BIE posed on the infinite boundaries of the domain. The perfectly matched layer (PML) coordinate stretching that is used to render propagating waves exponentially decaying, allows for the effective truncation and discretization of the BIE unbounded domain. We show through a variety of numerical examples that, despite the logarithmic growth of the complex-scaled Laplace's free-space Green's function, the truncation errors are exponentially small with respect to the truncation length. Our formulation uses only simple function evaluations (e.g. complex logarithms and square roots), hence avoiding the need to compute the involved water-wave Green's function. Finally, we show that the proposed approach can also be used to find complex resonances through a \emph{linear} eigenvalue problem since the Green's function is frequency-independent

Diffraction par un obstacle situé dans un réseau de plaques semi-infinies

International audience ; L'objet de cette Note est de proposer une méthode pour l'étude de la diffraction par un réseau de plaques horizontales semi-infinies localement perturbé par un obstacle. La méthode proposée couple une équation variationnelle posée dans un domaine borné entourant l'obstacle et une équation pseudo-différentielle écrite sur la droite verticale située à l'extrémité des plaques. Après avoir donné une formulation variationnelle du problème, on montre que celui-ci relève de l'alternative de Fredholm, en dehors des fréquences de résonance du réseau

Mesh requirements for the finite element approximation of problems with sign-changing coefficients

International audienceTransmission problems with sign-changing coefficients occur in electromagnetic theory in the presence of negative materials surrounded by classical materials. For general geometries, establishing Fredholmness of these transmission problems is well-understood thanks to the T-coercivity approach. Moreover, for a plane interface, there exist meshing rules that guarantee an optimal convergence rate for the finite element approximation. We propose here a new treatment at the corners of the interface which allows to design meshing rules for an arbitrary polygonal interface and then recover standard error estimates. This treatment relies on the use of simple geometrical transforms to define the meshes. Numerical results illustrate the importance of this new design